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Abstract. The thermodynamics of finite open antiferromagnetic XXZ chains 
is studied using field theory, Bethe Ansatz and quantum Monte Carlo methods. 
For the susceptibility a parameter-free result as a function of the number of sites 
L and temperature T beyond the scaling limit is derived. The limiting cases 
T/J ^ 1/L (J being the exchange constant), where the boundary correction 
shows a logarithmically suppressed Curie behaviour, and T/J <C 1/L, where a 
crossover to the ground state behaviour takes place, are discussed in detail. Based 
on this analysis we present a simple formula for the averaged susceptibility of a 
spin-1/2 chain doped with non-magnetic impurities. We show that the effective 
Curie constant has a highly non-trivial temperature dependence and shows scaling 
in the low-temperature limit. Finally, corrections due to intra- and interchain 
couplings and implications for experiments on Sr2Cui_3;Pda;03+a and related 
compounds are discussed. 
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1. Introduction 

The spin-1/2 Heisenberg chain is probably the first and most studied example 
of a strongly correlated quantum system. The first foundations towards exact 
solutions date back to the early days of quantum mechanics [1] and large scale 
numerical simulations were already performed in the early sixties [2]. Due to the 
interesting behaviour of this simple model and the connection with many low- 
dimensional antifcrromagnctic materials, the spin-1/2 chain is still heavily studied 
today. Some of the most fundamental properties, such as the exact amplitude of 
the power-law correlation functions [3-5] and the full temperature behaviour of the 
susceptibility [4,6], have only been established after the discovery of high temperature 
superconductivity has lead to renewed interest in low-dimensional antiferromagnetism. 
More recently, corrections from impurities and boundaries have come into focus [7-26] , 
which will also be the topic of this paper. 

Boundary thermodynamics are especially relevant for experiments [27-31] since 
even the most carefully prepared samples contain imperfections which effectively cut 
the chains, and corrections to the thermodynamic limit become especially large at low 
temperatures. In order to analyse the pure and the impurity contributions it is useful 
to define the finite size corrections 

5Ffs{L) = F{L) - LA„,, and 5xfs{L) = x{L) - iXb„ik (1.1) 

where F{L) {x{L)) is the total free energy (susceptibility) of a system with L sites and 
/bulk (Xbuik) is the free energy (susceptibility) per site of an infinite pure system. In the 
limit of L ^ oo, SFps and Sxfs become L-independent thermodynamic boundary 
contributions, which will be denoted by Fg and Xb and discussed in detail later on. 
However, in general L will also take on small values in experimental systems according 
to a random distribution. 

It is well understood that finite chains with an odd number of sites L have doublet 
ground states, which leads to a diverging Curie susceptibility at temperatures below 
the first excitation energy T/J ^ 1/L [7]. Recently, it has been established that the 
boundary susceptibility Xb{T) = Sxfs{T, L oo) is also divergent in the limit T ^ 0, 
albeit with a Curie factor which has a logarithmic temperature dependence [10, 11]. 
For experimental systems it is desirable to know the complete crossover between 
the two limits, which has to be averaged according to a random distribution. If 
the impurity concentration is p and we assume a Poisson distribution, i.e., that the 
impurity positions are uncorrelated, then the average susceptibility becomes [12, 26] 

oo 

Xp =/ 1^(1 -P)'^x(i) = (1 -p)Xbuik +P <5XFS, (1.2) 

where Sxfs = pX^l(1 — p)^Sxfs{L) and we have assumed that each impurity 
effectively removes one site. The typical approximation in experimental papers to 
estimate the impurity concentration by the low-temperature Curie tail corresponds 
to setting Sxfs = 1/(87") (assuming that each chain segment with odd length acts 
effectively like a free spin-1/2). This, however, docs not capture the full complexity 
of the problem and might, in particular, lead to an underestimate of the impurity 
concentration by an order of magnitude. 

Several experiments have been studying the prototypical spin-1/2 chain 
compound Sr2Cu03 [27,28,30,31] and the doped system Sr2Cui_a;Pdx03 [29] where 
palladium serves as a non-magnetic impurity. The analysis of the susceptibility even 
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in the undoped system has been hampered by large Curie-hke contributions at low 
temperatures making a detailed test of the theoretical predictions for the susceptibility 
of the pure system [6] difficult. In [27, 28] it has been shown that these contributions 
can be largely suppressed by annealing the samples. It therefore has been suggested 
that the Curie-like contribution is caused by excess oxygen. Each excess oxygen atom 
would then dope two holes into the chain. A detailed theoretical analysis of the 
susceptibility data [26] has shown that it seems to be possible to consistently describe 
the data by assuming that these holes are basically immobile and therefore effectively 
cut the chain into finite segments. A more detailed understanding of these experiments 
is desirable because it might also shed some light on the general issue of oxygen doping 
in cuprates, and in particular, in the high-temperature superconductors. 

Our paper is organised as follows: In section 2 the effective field theory description 
of the spin-1/2 chain is reviewed. In section 3 we calculate the free energy for an open 
anisotropic Heisenberg chain of finite length at finite temperature by field theory 
methods beyond the scaling limit. Next, we consider the limit L ^ oo with T/J 
fixed in section 4 where 5Fps and 5xfs as defined in (1.1) become well defined 
boundary quantities. In section 5 we then discuss separately the isotropic model 
which is the most interesting case from an experimental point of view. In the limit, 
T/J <C l/L, considered in section 6, a crossover to the ground state limit takes 
place and logarithmic corrections to the scaling limit result are studied in detail. In 
section 7 we then derive a simple formula which can be used to fit the experimentally 
measured averaged susceptibility Xp- We also discuss here the temperature dependence 
of the effective Curie constant C ^ TSxfs- In any real material there will be 
residual couplings bridging between the chain segments (intrachain couplings) as well 
as interchain couplings which will usually lead to some sort of magnetic order at 
low temperatures. The effect of such couplings will be discussed in section 8. A 
comparison between our theoretical predictions and experiments on Sr2Cui_xPdx03 
[29] is presented in section 9. Finally, we conclude and give a short summary in section 
10. 



2. Effective field theory for the spin-1/2 chain 

The Hamiltonian for the spin-1/2 Heisenberg chain with open boundaries in an 
external magnetic field h is given by, 

L-l L 

H = JY,[SiSt^, + Sl^l, + A5f -hY^St. (2.1) 

i=l i=l 

with J > for an antiferromagnetic system. Here, we have introduced an 
exchange anisotropy A which is convenient to obtain a low-energy effective theory 
by bosonisation. Such a description is by now well established [4, 32] and numerical 
simulations have shown its applicability for i > 10 and low temperatures T < 0.1 J, 
which is also the most interesting regime for experiments. In the massless case, 
— 1 < A < 1, the low energy fixed point of (2.1) is the Tomonaga-Luttinger liquid, 
which belongs to the universality class of the Gaussian theory with the central charge 
c= 1. 

Ignoring irrelevant operators the Hamiltonian is then equivalent to a free boson 
model 

Hq=^- j\x [tf + {d^(j)f] - h^^'^dx d^cp, (2.2) 
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where a' is the lattice constant. </> is a bosonic field and 11 = v~^dt(t> the conjugate 

momentum obeying the standard commutation rule [0(x),n(x')] = i5{x — x'). The 
dependence of the Luttinger parameter K on anisotropy A is known exactly from 
Bethe Ansatz 

K = [l--cos-^{A)]-^, (2.3) 

TT 

as well as the velocity [33] 

In the following we will set the lattice constant a' = 1. Factors of a' can be easily 
restored in the final results by dimensional analysis. Using the mode expansion for 
open boundary conditions (OBCs) 

n— 1 ^ 

the Hamiltonian can also be expressed in terms of bosonic creation and annihilation 
operators 

CO 

H, = + ^ ^n(aia„ + 1/2) - hS,. (2.6) 

n=l 



The zero mode operator Sz = y Jq dx dx<j) measures the total magnetisation of 
the chain and is therefore quantised to take on integer values for chains with an even 
number of sites L and half-integer values for odd L. 

Up to this point the effective model is easily solvable exactly. The partition 
function can be obtained directly by summing over all eigenvalues 

oo 

Z, = Y^e-^^'^^'^ n (2-7) 



e ( e~T^,-i^ 1 I I |2sinh 



OO 

" irv 



n[2 



2Ty J-J- L \2LT 

n=l 



where 0{q,u) is the elliptic theta function of the third kind = 9^{q,u) = 
Y^'^^-oo'i^ g»27iM fQj. integer Sz (even L) and of the second kind 6 = ^2(5, w) = 
Er=-oo g("+V2)'e^(2n+i)u for half-integer Sz (odd L). 
The free energy per site is then given by 

A = _|{,.,(.-*,_4)_|;to[2sinh(^„)]} . (2.8) 

ln6i(e-^,u) (2.9) 



and the susceptibility by 

52 „ 1 Esz^ze-^^' 1 9 



u=0 

Note, that xo has a simple scaling form as a function of LT [7] . The following limiting 
cases are of interest 

-Ur ^^P [-^h] ^'^h ^ 0, i even 

Xo= { ik LT/v ^ 0, i odd (2.10) 

± exp {-^Kif) LT/v - 00 . 
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Figure 1. The scaling behaviour of the susceptibiUty as a function of LT/v in 
Eq. (2.9) for the isotropic case A = 1, K = 1, v = 7r/2 (black solid lines) for L 
even (upper panel) and L odd (lower panel), respectively. The symbols denote 
quantum Monte Carlo data for chains of different lengths, obtained within the 
SSE loop algorithm framework [34]. 



Here the plus (minus) sign in front of the second term in the last line corresponds to 
L odd (even), respectively. This term represents the leading finite size correction to 
the thermodynamic limit result xo — K/(2'kv). In the limit T/v ^ 1/L on the other 
hand, the susceptibility for a chain of fixed even length vanishes exponentially because 
it locks into its singlet ground state, whereas xo for a chain of odd length exhibits a 
1/T divergence due to its doublet ground state. 

If we were only interested in the qualitative behaviour we could stop at Eq. (2.9), 
which in some cases is sufficient for agreement with experimental data [8]. However, 
the comparison with quantum Monte Carlo data in Fig. 1 clearly shows sizeable 
deviations from the scaling limit result. For more accurate predictions and the 
understanding of impurity corrections it is therefore necessary to take the leading 
irrelevant operator into account, which is given by [4, 32] 

SH = X [ dx cos{V8nK(t>). (2.11) 

The operator cos{^/ Sn K (f)) has scaling dimension 2K so that this perturbation becomes 
marginally irrelevant at the isotropic point K = 1. The next-leading irrelevant 
operators at zero magnetic field are of the form ~ {dx4>)'^ and will be neglected in 
the following unless it is explicitly stated otherwise. The coupling strength A has been 
obtained exactly by Bethe ansatz [4] and is given in a particular normalisation as 
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discussed in detail later on by 

KV(K) simv/K 



A 



7rr(2 - K\ 



r 1 



2K-2 



2V^r 1 



K 



2K~2 



2K-2 



(2.12) 



At the isotropic point A = 1 {K = 1) where the interaction (2.11) becomes marginal 
the amplitude A as given by (2.12) vanishes. Here A has to be replaced by an RG- 
improved coupling constant, and, in general, becomes a function of the different scales 
T, L, and h. This will be discussed in detail in section 5. 

Generally, in systems with boundaries, there may be boundary operators in 
addition to bulk interactions. However, as was pointed out in Ref. [13], in the absence 
of symmetry-breaking external fields at boundaries, there are no relevant boundary 
operators for the anisotropic Heisenberg chain. The leading irrelevant boundary 
operator has scaHng dimension 2 | and is of the form 6{x){dx(j))'^ [7]. This operator 
effectively leads to a replacement of the length L of the system by some effective length 
L' as will be discussed in the next section. 



3. Thermodynamics for a finite chain beyond the scaling limit 

To calculate the free energy and susceptibility for an open XX Z chain with length 
L at finite temperature T beyond the scaling limit (2.8,2.9), the leading correction 
(2.11) has to be taken into account. As this operator is irrelevant for < A < 1 it is 
sufficient to use perturbation theory. This anisotropic regime will be considered here. 
In first order perturbation theory wc obtain 

Sfi = ^J^ dx (cos (y8^(}>) J> . (3.1) 

This term has been considered in Refs. [10, 11, 25] in the limit L — > oo and gives the 
leading contribution to the boundary free energy. A bulk contribution, i.e. a term 
which scales linearly with L at large L, occurs in second order perturbation theory in 
A [4] (see also Appendix A). 

The expectation value can be split into an S'x-part (zero mode) and an oscillator 

part 



exp (±iV8TTK(P^ ^ = ^exp (±W8TTK(Pj ^ ^exp (±W8wK(j)^ ' 

exp(±iV8^<l>))^ exp(-47rii:(.^</.)o3c.) , (3.2) 



where we used the cumulant theorem for bosonic modes in the second line. Using the 
mode expansion (2.5) we find 



e L e Ki'T e t 
— 

KLT Q T 

and 



(exp [±iV8^<t>)) = _^ ^ (3.3) 



1 Note that the critical dimension for a boundary operator in (1 + 1) dimensions is 1, whereas it is 

2 for bulk operators. 
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The zero temperature part of (3.4) is divergent and we have to introduce a cutoff a 
with dimensions of length, of order a lattice spacing. Doing so and using X^^^ /I = 
— ln(l — z) for l^l < 1 leads to 



^ — ' TTt ZTT 



(3.5) 



1=1 



+ 



i-ln[(l-e2^*^/^e-"-/^) (l - 



g-27ria;/Lg-a 



Writing the exponential factor of the finite temperature part in Eq. (3.4) as a geometric 
series we find 



1 n=l 



-j^ ^nixl/L _|_ ^—27rixl/L 



21 



-■Kvnl/{TL) 



^ oo 

= i^{_ln(l-e-'^™/(^^)) 



(3.6) 



i.n[(l- 



^2-Kix/L --Kvn/iTL) 



-2-!rix/L^-TTvn/{TL)^ 



Combining Eqs. (3.4, 3.5, 3.6) gives 

cxp{-47ri\:(</-(/))„„,.} (3.7) 

-^^^ ^1 — e^'^*'^/'^e~"'^/^j ^1 — e~^'^*^^''"e~"'^^^ j ^ 

' ^-^^ _ Q2mx/LQ-vvn/{TL)^ _ g-27rix/ig-7rt>n/(TL)y 

(l_e-W(Ti.))2 



n 



-2K 



2an 

IT 



2K 



&K l^-^v/{TLY 



X Sin 



2^ /7ra\ ^ 



V L / 6»P (2^,e-'^''/(2r^)) 

In the last step we have written the oscillator part in terms of the Dedekind eta- 
function 

oo 

n=l 

and the elliptic theta-function of the first kind 

oo 

^i(m,9) = 2gV4 sin m[| (1-292™ COS 2u + /") (l-g2") . (3.9) 

n=l 

Prom (3.3) and (3.7) we now directly obtain the correction to the scaling form of the 
free energy (2.8) 

6f, = 2x(-) v'^'ie-^) dy-^ 

\LJ 'Jo ['iry,e ^tl) 

-I -K 



-K 



X 22^sin2^ (Try) [(l - e'^^^v ^-oi^ / l'^ ^ ^.c. 



(3.10) 
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g{y,q,w) = 55—5 (3.11) 

and A = a^^ X. It is actually A which is given by Eq. (2.12), but we will identify both 
in the following. Note, that the integral in (3.10) is convergent for K < 1/2. In this 
case the cutoff a can be dropped in the last line in (3.10) which then becomes equal 
to one. For K > 1/2 we can isolate the cutoff independent part by subtracting the 
Taylor expand of the integrand up to sufficient order in y, and then take a ^ 0. 

We will discuss this in detail in the following for the correction to the susceptibility 
which is readily obtained from (3.10) and reads 

'^Xi t2{l) ^ ' Jo ^ 01^ {7ry,e-^) 

X 2^^ sin^^ (Try) (^1 - e^'^'^'e""''/'^) x h.c. ^ (3.12) 
where we have defined a new function for the zero mode part 

.0 (.,.) = + j^-^. . (3.13) 

Because of the modified zero mode part, the integral is now convergent for X < 3/2 
and the last line can be set equal to one in this case again. For 3/2 < K < 5/2 we 
can obtain a convergent integral, i.e. the cutoff independent part, by subtracting just 
the first non-vanishing order in a Taylor expand in y and setting a ^ then. Noting 
that 

22^sin2^(7r2/)r?6^(e-f£) , , 

hm ; L = 1 3.14 

y^o (7ry,e-5Tr) 

the cutoff independent part of (3.12) for 3/2 < K < 5/2 is given by 



2A /^^2^^^; T .^70 .'y-^-'^ 2g2 ,,,,, 



where 



92 {q) = — ^ -^-2- • (3-16) 



Now we have to add again the first non- vanishing order in the Taylor expand in y but 
this time we keep the cutoff a giving us the non-universal contribution 

= '-^ Q'l (e-*) . (3.17) 

Shifting the upper boundary of integration to infinity, we can evaluate the integral 
and find 

'.r-^..(.-*)2f^„»-«. ,,18, 

The other non-universal correction stems from the irrelevant boundary operator 
~ 5{x){dx4>Y ■ Including this term into the Hamiltonian (2.2) is equivalent to replacing 
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the length by some effective length L ^ L' = L + a. Using this effective length in the 
exponentials in (2.9) and expanding to lowest order in a yields 

^X2 = ■j^f^92 (e-^) . (3.19) 

We see that this correction has the same form as (3.18). We therefore can consider a in 
(3.19) as an effective parameter incorporating both corrections. In the thermodynamic 
limit c/2 (e-"'"/^^^^^'^) K^T^L^ /{2Tr^v^) and 6x2 ^ Ka/{2TrvL). In this hmit 
we can compare the field theory result with a recent calculation of the boundary 
susceptibility based on the Bethe ansatz [24, 25] leading to 

a = 2^1/2 [ttK/{AK - 4)] / cos [7r/(4i^ - 4)] . (3.20) 

Eq. (3.12) for K < 3/2 or (3.15) for K > 3/2 taken together with (3.19, 3.20) is 
therefore a parameter-free result for the susceptibility of an open chain with length L 
at temperature T to first order in the Umklapp scattering. In [26] it has been shown 
by comparing with quantum Monte-Carlo data that this formula does describe the 
susceptibility of open chains for L > 10 and T/ J < 0.1 very well. Note, however, that 
the parameter a in (3.20) has poles at K = (4n + 3)/(4n + 2), n = 0, 1, 2, • • • with an 
accumulation point at K = 1. At these special points a contribution stemming from a 
different irrelevant operator will also be divergent while having the same dependence 
on temperature and length and both terms taken together will give something finite. 
In the limit L — » oo this has been investigated in detail in Ref. [25]. In particular, 
it has been shown in this limit that at K = 3/2 (A = 1/2) the two contributions 
(3.12) and (3.19) "conspire" to produce a logarithmic temperature dependence of the 
boundary susceptibility. 



4. Boundary contributions: The limit L — > oo with T/J fixed 

The boundary free energy can be obtained directly by starting with the integral (3.1) 
in the limit L ^ oo and using either boundary conformal field theory [10] or the 
thermodynamic limit result for the correlation function [11,25]. Here we want to show 
that our more general result (3.10) for a finite chain reduces to the known result in 
the thermodynamic limit. First, we consider the limit L ^ oo for the function g in 
(3.11). For a; ^ L we can replace the sums by integrals and find 

fif(^— ,e T^,eTj we cos I 1 cos I I (4-1) 

We see that g decays fast away from the boundary. If a; ^ L, on the other hand, we 
can set x — > i — .t and because Sz is integer or half-integer cos(47rS'2(iy — x)/L) = 
cos{A'K S zX / L) and we get exactly the same result again. In the thermodynamic limit 
the two boundaries therefore become independent and each of them yields the same 
contribution. To perform the thermodynamic limit for the oscillator part it is easiest 
to go back to Eqns. (3.4, 3.5, 3.6). We want to consider here only the case K < 3/2. 
In this case the cutoff in the second line of (3.5) can be dropped leading to 

CX).9/,;r\ -, ('r/'rv n /' 2tTx\ \ 1/2 ~1 
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Note that a will later be absorbed again into the coupling constant A — > A = a^^X. 

For the finite temperature part we can replace the momentum sum in (3.6) by an 
integral, allowing us to perform the sum over n ( "Matsubara sum" ) exactly 

^eeO-^ — — 



' 1 Jlikx I ^—2ikx 



— > In IH = — In — — — smh 



rt=l 

Thus we find in the thermodynamic limit 



(4.3) 



cos(^) 



Fb=L 6h{L ^ 00) = 2A / dx 2K ■ (4-4) 

h [^sinh(^)]^" 

This is the result obtained previously by boundary conformal field theory [10] and by 
using the result for the correlation function in the thermodynamic limit in (3.1) [11,25]. 
For ii" < 3/2 (A > 1/2) the integral is convergent and yields [25] 

FB = \MB{K + t^,l-2K)\(^-^Y (4.5) 

where B{x,y) = r(.x)r(y)/r(x + y). For K > 3/2 the integral (4.4) is divergent and 
a cutoff has to be introduced again. Then, the cutoff independent part of the integral 
is still given by (4.5) whereas the cutoff dependent terms are regular in h,T [25]. In 
the following discussion we will neglect these regular terms. Prom (4.4) the boundary 
spin susceptibility is easily obtained 

with il;'{x) = d\lj{x)/dx. Note that for 1 < < 3/2 (1/2 < A < 1), the boundary spin 
susceptibility xb shows a divergent behaviour ^ l/T'^^^^, as temperature decreases. 
This anomalous temperature dependence is also observed in the boundary part of the 
specific heat coefficient given by, 



-j {2K-l){2K-2)B{K,l-2K)^—^ . (4.7) 

We would like to stress that in the formulas (4.6) and (4.7) there is no free parameter, 
and the pre-factors are exactly obtained. These divergent behaviours are physically 
understood as follows. In contrast to the bulk Heisenberg chain in which the ground 
state is a spin singlet state, spin singlet formation in the vicinity of boundaries is 
strongly disturbed. For the susceptibility, for example, we can write xb = Xb i^) 
where Xb(^) follows from (4.4). Note, that the local boundary susceptibility has 
also an alternating part [18], which, however, does not contribute to xb- A detailed 
comparison of numerical data for the local boundary susceptibility (consisting of the 
uniform and the alternating part) and field theory has been presented in Ref. [24]. 
In Fig. 2, x^i^') is plotted for various temperatures showing that the spins near the 
boundary are more susceptible than the spins in the bulk. It should, however, be 
emphasised that the singular behaviours are not due to the presence of boundary 




Figure 2. The uniform part of the local boundary susceptibility Xb(^) following 
from Eq. (4.4) for K = 1.2 at different temperatures. The boundary susceptibility 
Xb is obtained by integrating over Xs(^) leading to a 1/J'3-2K divergence if 
K < 3/2. 



operators, but interpreted as a consequence of finite-temperature corrections caused 
by bulk irrelevant interactions. 

At zero temperature with a finite magnetic field, a similar singular behaviour 
appears in the field dependence of the boundary spin susceptibility given by, 



Xb(T = 0) = A — sin(7rX)r(3-2X) . (4.8) 



V 



The zero temperature susceptibility can also be derived from the Bethe ansatz exact 
solution by using the Wiener-Hopf method. This has been done in Ref. [24, 25] and 
agrees with the expression above. 



5. The isotropic case A = 1 

At the isotropic point Umklapp scattering becomes marginal. Instead of using (2.11) 
one can write the marginal perturbation in a manifestly SU(2) invariant form using 
the SU(2) current operators. The perturbation in this formulation becomes 

SH = J g[5||J,J, + ^(J+J_ + J_J+)]. (5.1) 

The small running coupling constants g\\ , g± fulfil a set of known Renormalisation 
Group (RG) equations [3,4] and are in general functions of the different scales, T,L 
and /i, involved. The perturbation (5.1) can be split into two parts: JzJz oc lP — {dx4>)^ 
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and J+J_ + J_ J_)- cx cos(v87r(/)). The first term can be absorbed in the free electron 
Hamihonian (2.2) by a rcscaling of the fields 11 and </). This leads to a rcnormalisation 
of the Luttinger parameter K l+g\\/2 + 0{g^) (see Appendix B for details). With 
this replacement in (2.9) we see that the susceptibility at the isotropic point in the 
thermodynamic limit is now given by 

X.u..(T) = l (^1 + ^ + 0{g^)^ , (5.2) 

a result first derived in Ref . [6] . The second term leads in lowest order to a replacement 
A — > — p_l/4. More generally, the results for the anisotropic model derived in the 
previous sections have to be expanded in powers of 1 — 1/K and then re-expressed in 
terms of the coupling constants g|| and g±^ [4]. For details the reader is referred again 
to Appendix B. Here we only want to give the main results. 

First, we consider the isotropic limit of our first order result for the boundary 
susceptibility (4.6). This leads to 

(1) _ g± 9\\9± ( 1 «'"(!) 



~ 12T 8T U 7r2 J ■ ^^'^^ 
The second order result derived in Appendix A also yields a contribution in quadratic 
order in g 

Xi'^«f|-0.11. (5.4) 

Here the factor 0.11 stems from a numerical evaluation of the integral in (A. 9) for 
K = 1. With g\\ g, g± ^ —g we therefore obtain up to quadratic order in g 

i2r 8r 



XB = ^,+-^ + ^[ 0.66- ^ ) +0{g^). (5.5) 



Here g = g{T) and is given by [4] 

l/g + ln{g)/2 = ln{To/T) , (5.6) 

with To = ^/n/2e^/^+'1 where 7 is the Euler constant. Note, that the scale in the 
logarithm, Tq, is chosen according to the expansion of the bulk susceptibility. Tq is 
non-universal and can be fixed differently to achieve the best possible agreement with 
the expansion of the boundary susceptibility at the isotropic point. However, this is 
not important for XB as it only influences sub- leading terms so we keep the scale as 
set by the bulk part in order to have a single g in the formula for x(L, T) discussed 
in the following. A more detailed discussion of this point is presented in Appendix B. 

In addition, we have added a contribution ^ Ka/{2T:v) to (5.5) with some 
constant a stemming from the boundary operator as in the anisotropic case. However, 
the constant (3.20) obtained by Bethe Ansatz in the anisotropic case does have an 
accumulation point of singularities at K = I making it impossible to extract the 
value at the isotropic point from this formula. If one considers, on the other hand, 
the integral for the boundary susceptibility derived directly by Bethe Ansatz for 
the isotropic case (see Eq. (36) in [24]), this task is difficult as well because the 
logarithmic contributions coming from the cut along the imaginary axis completely 
dominate. Indeed, at low-temperatures the constant contribution is less important 
than any of the logarithmic terms coming from the expansion of xb outlined above. 
Nevertheless, this constant becomes important again at higher temperatures (but still 
T/J<^ 1) as in the anisotropic case. We will therefore use it as a fitting parameter 
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Figure 3. Comparison between QMC results (dots) and tiie field theory result 
(5.5). The black solid line corresponds to taking into account the first order in g 
term only, the blue dashed line is the result with the second order term included, 
and the red dot-dashed lino is obtained by adding a constant Kajl-KV with a = 1.5 
[K = 1, V = 7r/2). The QMC data for XB have been obtained by subtracting the 
bulk susceptibility known from Bethe ansatz [6] from numerical data for xi^)- 
The chain length L has always been chosen such that TL/v > 4. In addition, the 
susceptibility has been calculated for lengths L and L + 1 and then averaged to 
further reduce finite size effects (see also Eq. (5.8) and Fig. 4 below). 



when comparing to the numerical data in Fig. 3. Note, that a now effectively also 
partly incorporates higher order logarithmic corrections to (5.5) as well as the constant 
coming from the boundary operator. Furthermore, we must in principle replace 
Ka/{2-Kv) — > (a/7r^)(l + g/2 + + • • •) (see Appendix B). However, because the 
logs dominate at low temperatures this merely leads to a small rescaling of a which 
we will ignore because a is an "effective constant" in any case. Solving Eq. (5.6) we 
can also write the boundary susceptibility as 

^ 1 A _ inin(ro/r) 

12T\n{To/T) \ 21n(To/r) 

This means that even in the limit T/v ^ 1/L we find a Curie-like contribution albeit 
with a "Curie constant" which depends logarithmically on temperature. 

If T/v becomes of order 1/L the susceptibility can no longer be split into a 
bulk and a boundary contribution. However, we can still subtract the known bulk 
contribution and define a quantity Sxfs as in Eq. (1.1). In a first approximation, 
we can use the thermodynamic limit result for xb, Eq. (5.5), and add the leading 
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finite size corrections stemming from tlie scaling part (second term in the last line of 
Eq. (2.10)). This leads to 

or 2^1 

SxFsiL, T) « xb{T) ± ^(1 + g) cxp[-2Lr(l + g/2)] (5.8) 

where the plus (minus) sign refers to chains of odd (even) length and we have replaced 
K ^ 1 + g/2 with g determined by Eq. (5.6). A comparison between this formula 
and QMC data is shown in Fig. 4. Note, that Eq. (5.8) implies an almost perfect 




100 200 300 400 500 600 

L 

Figure 4. QMC data for the susceptibility <5xfs s-*- temperature T = J/30 as a 
function of chain length. These data have been obtained by subtracting the bulk 
susceptibility known from Bethe ansatz [6] from numerical data for x(-L). The red 
squares (blue circles) correspond to even (odd) chain length, respectively. The 
lines represent the approximation (5.8). 

cancelation of finite size corrections when calculating the susceptibility for lengths L 
and L + 1 and then taking the average. This property has been used to improve the 
numerical data for xb presented in Fig. 3. 

To obtain the result for x(I/,T) at the isotropic point in general, we have to 
consider both parts of the perturbation (5.1): First, we have to replace K ^ 1 -\- g/2 
in the exponentials of (2.9). For the correction (3.12) we can only obtain a result to first 
order in g because we do not know the analytical solution of the integral and therefore 
cannot easily expand it around K = 1. The result to first order in g is readily obtained 
by just replacing A ^ g/4 and evaluating the integral for K = 1. Now, however, the 
running coupling constant g depends on two scales: the length L and temperature T. 
There is no general solution of the renormalisation group equations if two different 
scales are involved. At low enough energies, however, the smaller length scale will 
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always dominate and the running coupling constant becomes [4] 

l/g + lnig)/2^lTi(^^/2/^e^/'^+''mm[L,v/T]j . (5.9) 

In addition, the constant a stemming from the boundary operator is still present and 
we will use it again as a fitting parameter. Note, that a now will be different from the 
one obtained by fitting the numerical data for xb as it also partly incorporates the 
logarithmic terms now taken into account only to lowest order. 

To summarise, the susceptibility for the isotropic case is given by x(-^> 2^) = 
Xo + <5xi with the replacements X — > 1 + g{L, T) /2 + a/L in the exponentials of (2.9) 
and A g{L,T)/A in (3.12). A comparison between this formula and QMC data for 
chains of different lengths is shown in Fig. 5. Here we find that a « 4 works best 




T 

Figure 5. Comparison for A = 1 between QMC results (symbols) and the field 
theory result (lines) with a = 4 for chains of different length L. 

describing the numerical data over a wide temperature range and for different even 
and odd chain lengths very accurately. § 

The disadvantage of this formula for x{L,T) is, that the integral in (3.12) has 
to be evaluated numerically for each length and temperature considered. In the next 
section we will derive a much simpler formula in the limit T/J <^ 1/L. We find 
empirically that this formula describes the susceptibility over a much larger range of 
temperatures than anticipated making it useful to fit experimental data without the 
necessity to evaluate the integral in (3.12) and without having a free fit parameter. 

§ In [26] we used a = 5.8. This difference is due to the fact that in this article we expanded in 
a, i.e., we used (3.19) as well as a second order solution for (5.6). Here we put a directly into the 
exponentials of (2.9) and solve (5.6) numerically. 
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6. The ground state limit: LT/v^O 

Next, we want to consider the limit LT/v 0. We have already shown in Eq. (2.10) 
that the scaling part of the susceptibility shows a Curie-like divergence if the chain 
length is odd and an exponentially suppressed behaviour if it is even. For the 
first order correction (3.12) we will concentrate again on the case K < 3/2 where 
the cutoff can be dropped. In the limit LT/v — > we find r] (e~ tt ) — > e~ ^irr ^ 
6i (tt?/, e~2¥r) — » 2e~^n: sin(7rj/) and 



9o (y, e 



L odd 

— 2cos(47rj/)e~'!^ + 2e~'!f^ L even 



(6.1) 



Here we have expanded go only to first power in the small parameter e~ t^tt . For L 
odd the first non-vanishing contribution is second order. As the scaling part shows a 
power-law divergence in this case, any exponentially small corrections can be safely 
neglected in any case. For L even, on the other hand, they are important as the scaling 
part is also exponentially small. In the even case we can now evaluate the integral 
and find 

L odd 



Sxi{LT/v^O) 



2A /ttN 

T2 II; 



2-fS' r(3-2_R-) 

r(2-_ff)r(3-i<:)* 



Finally, lets consider the contribution 6x2- For L odd 52(6 klt") 
whereas 52 (e""^^) — > 2e~T^ for L even. This leads to 



^^^ = { 2.. 
For the susceptibility in the limit LT/v 



KTL ■ a 



Lodd 

L even 
we therefore find 



(6.2) 

- L even 
in lowest order 

(6.3) 



X(LT^0) = 



1 

iLT 



1 + A 



L2K-1T 



L odd 



L even 



(6.4) 



where /? = r(3 - 2iC)/[r(2 - K)V{2> - K)]. So corrections to the scaling limit result 
in first order in e~Tr are only present in the case of even chain length. However, 
even in this case these corrections are suppressed by additional powers of 1/{LT). 
So one might think that Eq. (6.4) is of purely academic interest. However, as we 
will show below, it enables us to derive a simple formula for the susceptibility of the 
isotropic chain which works over a much larger range of temperatures and lengths 
than anticipated. 

In the isotropic case we have to replace K ^ l-|-5i||(L)/2 in the scaling part (2.9). 
Furthermore A —g_\_{L)/4:, so that (6.4) at the isotropic point reads 

L odd 

(6.5) 



1 

ALT 



X{LT ^ 0) 



L even 



where v = 7r/2 has been employed. Note that with g\^ g, g^^ —g the replacement 

K —i l + g||(L)/2 in the scaling part and the first order contribution in A ^ —g±{L)/4: 

both produce exactly the same correction j^g{L) in the brackets in (6.5). Here 
g{L) is given by (5.6) with L being the relevant scale in this limit. We want to 
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stress that both terms yielding exactly the same contribution in the Umit considered 

here contribute very differently in the thermodynamic Umit. Whereas the gp-part 
leads to a correction to the bulk susceptibility, the gx-part determines the boundary 
susceptibility. 

Ignoring the contribution ^ a which comes from the boundary operator and is 
suppressed by an additional power of L we can therefore obtain exactly the same 
expansion by replacing 1/K ^ 1 — g{L) in the exponentials of (2.9). I.e., we can 
absorb the first order correction 5xi into the scaling form. This leads to the formula 



1 



5f exp 



-^{l-g{L))Sl 



This formula should be valid in the; limit LT /v — > when only the ground state and the 
lowest excited states contribute to the sum over Sz- One easily verifies that Eq. (6.5) 
is indeed reproduced in this limit: For L odd we take only the states Sz = ±1/2 into 
account and the correction 1 — g{L) cancels out. For L even we consider the ground 
state Sz = Q and the first excited state Sz = ±1- Expanding in g{L) then leads to 
(6.5). 

So starting from our general result for the susceptibility of finite open chains we 
have found that the correction to the excitation energy of the lowest excited states 
with 5^ ± 1 in the isotropic case is given by i? — !■ {itv/L)[\ — g{L)\. This agrees with 
the findings in Rcf. [14]. Corrections to higher excited states are not of such simple 
form, and, in general, it is therefore not possible to absorb these corrections into the 
scaling part. 

Empirically, however, we find that (6.6) works quite well over a much larger 
range of lengths and temperatures than one might anticipate. In Fig. 6 we compare 
this formula again to QMC data. Although the agreement at higher temperatures 
where the susceptibility becomes constant is not perfect, the errors in this region are 
only about 3%. 



7. The averaged susceptibility and the effective Curie constant 

In this section we want to discuss the averaged susceptibility of the isotropic Heisenberg 
chain for a concentration of chain breaks p as defined in Eq. (1.2). Here we have 
assumed that we have a completely random (Poisson) distribution of impurities. 
However, we could obtain Xp for any other distribution as well. Xp can be calculated 
by using our field-theory results in section 5 for L > 10 and T/J< 0.1 in combination 
with exact diagonalisation or QMC results for shorter chains. This method has been 
used in Ref. [26] to analyse susceptibility data for Sr2Cui_a;Pdx03+5 [29]. Based on 
our analysis of the limits LT/v and LT/v ^ oo in the previous sections we want 
to derive here a simple formula which allows to determine the concentration of chain 
breaks by fitting the measured susceptibility. 

We have seen that there are no corrections to the scaling limit result in the limit 
LT/v for odd chains. For even chains the susceptibility is exponentially small 
in this limit in any case so that even chain segments practically do not contribute to 
the average in (1.2). In the opposite limit we can write xi^) ~ -^Xbuik + Xb- We now 
assume that the crossover occurs at some length Lc = "yJ/T where 7 is a crossover 
parameter which we expect to be of order one. For the average susceptibility we can 
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Figure 6. Comparison for A = 1 between QMC results (dots) and the field 
theory result (6.6) in the limit LT/v — > (lines). 



therefore write 

2 Lc oo 



AT 

L odd 

4T2-P 



■pY 



(7.1) 



(l _ (1 _ p)7/T) + (1 _ py/T [(l - p + ^) + PXB 



and we use the field theory result for the bulk susceptibility [4] 



Xb 



fi(T) , 3p3(r) 



(7.2) 



2 32 

where the T^-term stems from the irrelevant operators with scaling dimension 4. The 
boundary susceptibility xb is given by (5.5) with a = 1.5 as in Fig. 3. A comparison 
between QMC data for Xp and (7.1) with 7 = 1 for various impurity concentrations 
is shown in Fig. 7. 

The impurity part and an effective Curie constant C can then be defined according 
to Eq. (1.2) by 



pSxFS = Xp - (1 - p)Xh 



C 



(7.3) 



This effective Curie constant is shown in Fig. 8. Note, that only in the extreme low- 

tcmperature limit the usual Curie constant C l/4-(l — p)/(2 — p) w 1/8 (half 
of the chains contribute a Curie constant of 1/4 for <C 1) is recovered. However, 
C(2 — p)/(l — p) as a function of T/p shows an almost perfect collapse (scaling) for 
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Figure 7. Txp for p = 0.01,0.03,0.05,0.07,0.1 from bottom to top. The black 
dots are obtained by using QMC data for x{^) ^^'^ the lines represent the fitting 
formula (7.1) with 7 = 1. 




Figure 8. The effective Curie constant as defined in Eq. (7.3) for p = 
0.01, 0.03, 0.05, 0.07, 0.1 (from bottom to top) as a function of T/p. 
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different impurity concentrations if T/p < 1/2. At higher temperatures this scahng no 
longer holds because the boundary contribution xb is not a scaling function in T/p. 
The non-trivial temperature dependence of the effective Curie constant, in particular 
the finite temperature minimum, and the scaling at temperatures T <p/2 provide a 
way to test this scenario experimentally. 



8. Inter- and Intrachain couplings 

In any real system there are residual couplings between the chain segments in 
the principal chain direction (intrachain couplings). In addition, there are also 
perpendicular couplings between the chains (interchain couplings) which at low enough 
temperatures will destroy the one dimensionality of the system and usually induce 
some sort of magnetic order. In Sr2Cu03, for example, the nearest-neighbour 
Heisenberg coupling along the chain direction (b-axis) is estimated to be J ~ 2200 K 
whereas the couplings along the other axis are orders of magnitude smaller (J^ 5 
K, Jc ^ 10~^ K) [28]. Furthermore, it has been pointed out that the next-nearest- 
neighbour coupling along the chain direction is probably not that small, J2 ^ 140 
K [35]. This irrelevant coupling will only cause small corrections to the susceptibility 
of an isolated chain segment by changing the marginally irrelevant coupling constant 
A. However, it induces a coupling between different chain segments even when the 
nearest-neighbour exchange J is absent due to a non-magnetic impurity. Note, that 
there is no frustration in this case and J2 just acts as an effective coupling J' between 
the chain ends as shown in Fig. 9. Although the next-nearest neighbour coupling 



J' non-magnetic impurity 




Li 11 L2 



Figure 9. Two chain segments with length Li and L2, respectively, coupled by 
a weak intrachain coupling J' ^ J. 



across an impurity can differ from the one across a magnetic ion, we will assume that 
it is still of the same order of magnitude, i.e., J' ^ J2. 

The susceptibility of the weakly coupled chains is given by 

X=^(^{p^St, + i:Sl}j ^ (8.1) 

with L = Li + L2- Here Sf^ (S'l^) denotes the z-component of the spin on chain 1 (2) 
at site i (j), respectively. For J' ^ J we can treat this coupling using perturbation 

theory. In zeroth order this lead to x'"^ = {Lixf'' -I- £2X2°^)/-^- In first order only the 
longitudinal part of the coupling contributes and we obtain 

^^-'1^2 ^E^H^n) (Es^sS^^I 

IT' IT' 
- -JJf2 {SL.St,) {SL.S^,) - -— X11X21 (8.2) 
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where Xu = < 'S'fi > /{dh) is the susceptibility of the boundary spin 
and Sf^^^ = Y^^Li^ii {I = 1,2), respectively. The extra factor of 1/T arises due 
to time-translational invariance and the imaginary time integral in the first order 
perturbation theory formula. Because of the broken spatial-translational invariance, 
the local susceptibility is position dependent and the boundary spins are more 
susceptible than the spins deep inside the chain. These effects have been studied 
in detail in Refs. [18,24]. As we are here just interested in an order of magnitude 
estimate of the effect of intrachain coupling on Xp we can replace the susceptibility 
of the boimdary spin xii by the average susceptibility per site xj'^''. This leads to 

dx ~ — 2(J'/L)xi'''x2°^- Assuming further that for an impurity concentration p all 
chain segments have length L 1/p we find 

SXp ~ J'PXI ■ (8.3) 
At T ^ we have Xp ^ P/{^T). The correction (8.3) therefore becomes of order 
Xp at T ^ >^'P^-II Foi" temperatures T J'p^ intrachain coupling can therefore be 
neglected. 

Although first order perturbation theory can give us the temperature scale 

where intrachain coupling becomes important, it is not sufficient once this scale is 
approached. This can be seen by a scaling analysis. We have ~ dx4> with scaling 
dimension one. In n-th order perturbation theory in the intrachain coupling J', the 
longitudinal part of the coupling therefore yields x*-"'' ~ {TL)~^{J' /TL"^)^ . We can 
therefore write the susceptibility in scaling form 

Replacing again L 1/p we find for the averaged susceptibility 



This means that at temperatures where intrachain coupling J' can no longer be 
neglected, T ~ J'p^, perturbation theory in J' breaks down. At temperatures 
T <^ J'p^ < Jp, on the other hand, we can obtain an effective model by replacing 
all segments of odd length by 5 = 1/2 spins. The couplings between these effective 
S = 1/2 will then be random with a very broad distribution. According to Fisher [36] 
we might therefore expect that a system of weakly coupled chain segments renormalises 
at very low temperatures to the infinite random fixed point. In any realistic system, 
however, there will be also interchain couplings between the chain segments preventing 
the system from reaching this fixed point. 

An interchain coupling between two chain segments is shown in Fig. 10. In first 
order perturbation theory we obtain in this case 



<5x(^^ = - 7^ V ( V Sl.Sl, ) I V Sl,SS^,_, ) (8.6) 




S2,n+l-k) — X] Xl,fcX2,n+l- 



II Note that there is a factor p missing in Ref. [26]. 
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Figure 10. Two chain segments with length Li and L2, respectively, have an 
overlap of n sites and are coupled by an interchain coupling J. 



where xi,k denotes the local susceptibility of the /-th chain at site k. To simplify 
matters, we consider the case of two chains with equal lengths L\ = L2 = L which 
completely overlap, i.e., n = L. In this case Eq. (8.6) becomes 



8x^^) = - 



L ^ 

fe=i 



Xl,feX2,fe 



(8.7) 



At T <C we can again set xi.k ^ Xp ^ p/T which leads to SxjP ^ Jp^ /T^ . This 
becomes of order Xp at temperatures T ~ Jp. As J is of order of the Neel temperature 
Tjv this scale is completely irrelevant. However, for a system with open boundaries the 
local susceptibility has also a staggered part [12,18] which will contribute to (8.7). This 
is a consequence of the S^-operator having a staggered part, ~ c(— 1)*^ cos •\/27r0, 

where c is a constant. For a chain with odd length in the ground state limit, T <C 1/i, 
we find [12] 



± 



qz 

tot'-'k,st. 



±) « (-1)* 



T V 2L 



nk 



(8.8) 



where |±) is the ground state with S^g^ = ±1/2, respectively. This leads to 
^X^^^ ^ J /{LT"^) Jp/T"^ and therefore dominates compared to the contribution 
originating from the uniform part of the local susceptibility. It becomes of order Xp 
at temperatures T ^ J. 

Finally, we want to consider the "clean limit", T pv. In this case, the Neel 
temperature Tjv can be determined by treating the interchain coupling in mean field 
theory. This leads to the well known condition 

zJxstiTN) = 1 (8.9) 

where Xst is the staggered susceptibility, i.e., the staggered response to the effective 
staggered field originating from the other chains and z is the coordination number. 
Because Xst {T) ~ 1/T we obtain that Tjv ~ J. The same condition for the breakdown 
of one- dimensionality in this limit is obtained by considering the contribution of the 
staggered part of the 5^-operator to the susceptibility in second order perturbation 
theory in the interchain coupling. 

We therefore conclude that our results for the averaged susceptibility Xp the 
previous sections are applicable as long asT ^ Min(J,p^ J'). 



9. Experimental situation 



The best known realisation of the spin-1/2 Heisenberg chain is Sr2Cu03. Here copper 
is in a 3d® configuration and has spin 1/2. The Heisenberg couplings between these 
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spins are spatially very anisotropic as already mentioned in the previous section. For 
T 3> Ja 5 K interchain coupling can therefore be neglected and the magnetic 
properties are described by the one-dimensional isotropic Heisenberg model. 

The doped compound Sr2Cui_a:Pd2;03 has been studied by Kojima et al [29]. 
Here palladium has spin zero and therefore serves as a non-magnetic impurity. 
Impurity concentrations from x = 0.5% up to a; = 3% have been studied. This 
means that x^J2 with J2 140 K is small even compared to Ja- In this material 
the limit where our theory is applicable is therefore set by the interchain coupling, 
T > J„ « 5 K. 

The theoretical analysis of the susceptibility measurements on Sr2Cui_xPdx03 
is hampered by the fact that even the pure system Sr2Cu03 already shows Curie- 
like contributions at low temperatures [27, 28]. It has been observed that these 
contributions can be significantly reduced by annealing and two explanations have 
been offered: Excess oxygen might be present in the as grown compound. Most likely, 
each excess oxygen would then dope two holes into the chain leading to two Zhang- 
Rice singlet type states as in the high-Tc compounds. Assuming that these holes are 
relatively immobile, they effectively act as chain breaks. In this scenario the Curie- 
type contribution in the "pure compound" is caused by the mechanism described 
in this paper. Indeed, it has been shown [26] that the susceptibility data for the 
as grown crystal sample in [28] can be well described by assuming a concentration 
of chain breaks, p = 0.6%. For the powder samples studied by Ami et al. [27] an 
alternative explanation was proposed by Hill et al. based on the observation that 
Sr2Cu03 decomposes into Sr2Cu(0H)e under exposure to air and water [37]. The 
copper spins-1/2 in Sr2Cu(OH)6 are only weakly coupled according to Hill et al. and 
can therefore induce a true Curie contribution as opposed to the more complicated 
behaviour expected in the chain break scenario. It is, however, hard to see how 
this could also explain the Curie-like contributions observed on single crystals by 
Motoyama et al. [28] . We therefore believe that the excess oxygen scenario is the most 
plausible explanation in this case. 

The Sr2Cui_xPda;03 samples investigated by Kojima et al. have not been 
annealed so that we expect chain breaks due to excess oxygen in addition to the once 
caused by palladium. Furthermore, it can certainly not be excluded that there are 
indeed additional true paramagnetic impurity contributions caused by decomposition 
or other mechanisms. This makes a thorough theoretical analysis very difficult. In [26] 
such an analysis has been performed under the assumptions that true paramagnetic 
Curie contributions can be neglected and that excess oxygen will lead to immobile 
holes acting as additional chain breaks. The measured susceptibility x'""'' then 
consists of a temperature independent part x' due to core diamagnetism and Van 
Vleck paramagnetism and Xp coming from the Heisenberg chains with impurity 
concentration p. Here p = x + 5 has to be considered as a fitting parameter consisting 
of the nominal palladium concentration x and the concentration of additional chain 
breaks 6 due to excess oxygen. A best fit of the experimental data was obtained 
in [26] with the two parameters p, x' as displayed in Table 1. For the sample with a 
nominal Pd concentration of 3% our best fields has yielded p = 2.4%, i.e., less than 
the nominal concentration. One might speculate that the Pd ions cluster or that part 
go in at interstitial positions instead of replacing Cu ions thus reducing the number 
of chain breaks created. 

Here we want to go one step further and compare the effective Curie constant 
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Table 1. Concentration x of Pd ions in experiment compared to impurity 
concentration p and constant contribution x' yielding the best theoretical fit. 
The first line corresponds to the "as grown" sample of Sr2Cu03 from Ref. [28]. 

X (Exp.) p (Theory) x' [emu/mol] 



0.0 0.006 -7.42 X lO'^ 

0.005 0.012 -7.7 x 10-5 

0.01 0.014 -7.5 X 10-5 

0.03 0.024 -7.5 x 10"^ 



extracted from the experimental data 

C^'-" = - [X"''" - X' - (1 - p)Xw.] (9.1) 
P 

with the theoretical prediction. We will use the same values for the concentration 
of chain breaks p and the constant contribution to the susceptibility x' as in Table 
1. For Xbuik we employ the field theoretical result (7.2). The result is shown in 
Fig. 11. Although the agreement is far from perfect we see that the Curie constant 
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Figure 11. The effective Curie constant (9.1) extracted from susceptibility data 
for Sr2Cui_a;Pd2:03 [29]. The nominal Pd concentration is x = 0.5% (black 
circles), x = 1% (red squares), and x = 3% (blue diamonds). The solid lines 
denote the theoretical result (7.3) for p = 1.2%, p = 1.4%, and p = 2.4% (from 
bottom to top). 



extracted from experiment does indeed show a non-trivial temperature dependence. 
In particular, there is a finite temperature minimum although at higher temperatures 
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than predicted theoretically. The maximum at low temperatures and subsequent drop 
of the Curie constant, on the other hand, is clearly caused by the interchain couplings 
which can no longer be neglected in this temperature regime. Given that the procedure 
to extract the Curie constant is quite sensitive to the values of J, p, and x' used the 
agreement is quite remarkable supporting that the Curie-like contribution is mainly 
caused by an effective concentration of chain breaks p and that possible additional 
paramagnetic Curie contributions are relatively small. 

For a more detailed comparison with theory, it would be desirable to repeat 
the experiment with annealed samples, carefully estimating any residual true 
paramagnetic Curie terms in the pure system, so that the nominal Pd concentration 
directly corresponds to the concentration of chain breaks created. 

10. Conclusions 

We have studied the effect of chain breaks (non-magnetic impurities) on the 
thermodynamic properties of spin-1/2 chains by field theory methods. Using first 
order perturbation theory in the leading irrelevant operator (Umklapp scattering) we 
have derived a parameter-free result for the thermodynamic properties of a finite length 
XXZ-da&m with open boundaries beyond the scaling limit. We have shown that this 
result reduces to the previously known expressions for boundary contributions in the 
limit of infinite chain length. To obtain results for the isotropic case, we expanded the 
results for the anisotropic model in terms of 1 — l/JsT and re-expressed this expansion in 
terms of a small running coupling constant g fulfilling a set of known RG equations. For 
boundary contributions like the boundary susceptibility xb this leads to multiplicative 
logarithmic corrections. Whereas results in first order in g can be easily obtained by 
just replacing the Umklapp scattering amplitude A .9/4, the consistent calculation of 
the correction of order g" generally requires n-th order perturbation in the Umklapp 
scattering. For the boundary susceptibility we presented the result up to order g"^ 
using second order perturbation theory in A. 

For the susceptibility of a finite length isotropic Heisenberg chain we were able to 
derive a simple formula by studying the ground-state limit LT/v — > 0. In this limit the 
first order correction in the Umklapp scattering can be absorbed into the scaling limit 
result. Surprisingly, we find by comparing with QMC data that the simple parameter- 
free formula for x(L,T) obtained this way does work very well over a much larger 
range of temperatures and lengths than anticipated. This simple formula is useful to 
analyse susceptibility data for finite length Heisenberg chains. 

We then discussed the average susceptibility Xp for impurity concentration p 
assuming a Poisson distribution. Using the results for x(i,T) in the limiting cases 
LT/v and LT/v ^ oo we derived an effective fitting formula for Xp which 
reproduces QMC results with high accuracy. This formula can be used to fit 
experimental data thus allowing to determine the (effective) impurity concentration. 

Taking into consideration that any real material will have interchain J as well as 
longer ranged intrachain couplings J' which can bridge over a chain break we discussed 
the limitations of the theoretical analysis for Xp presented here. Our main result is 
that such couplings can be ignored provided that T Min(J,p^ J'). Note, that J' is 
suppressed by a factor so that even relatively large couplings between the segments 
can be ignored provided that the impurity concentration is small. In most materials it 
will therefore be the loss of one-dimensionality due to J which marks the point where 
our theory no longer applies. 
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For Sr2Cu03 this scale is about 5 K whereas J ~ 2000 K making it an 

ideal candidate to study the thermodynamics of impurities in a Heisenberg chain 
experimentally over a large temperature range. In this work we concentrated on 
analysing susceptibility data by Kojima et al. [29] for Sr2Cui_a;Pda;03. We extracted 
the effective Curie constant and showed that it agrees qualitatively with the theoretical 
prediction. The impurity concentration, however, seems to be different from the 
nominal Pd concentration. A plausible explanation is the presence of excess oxygen 
leading to additional chain breaks, a scenario already discussed in [27,28]. It might 
be possible to avoid these additional chain breaks by annealing which seems to make 
a more detailed comparison with theory feasible in the future. 

Our results might also shed light on some aspects of the high- Tc cuprates physics. 
In YBa2Cu306+5 we have quasi-one-dimensional CuO chains in addition to the Cu02 
planes common to all high- Tc cuprates. By changing the oxygen content 5 the charge 
concentration in the planes is changed. At the same time, this leads to the removal of 
spins from the chain. If we assume a Zhang- Rice singlet type state with immobile holes 
this again leads to chain breaks in the same way as in Sr2Cu03+5. Indeed, Friedel- 
type spin-density oscillations have been observed in YBa2Cu306.5 by NMR [38]. Such 
oscillations are expected if spin chain segments with open boundaries are present [18]. 
A thorough analysis of these data and possible new experiments on Sr2Cui_xPdx03 
or similar compounds might give us a better understanding of oxygen doping and the 
formation of Zhang-Rice singlet states in the cuprates. 
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Appendix A. Second order perturbation theory 

The free energy in second order in the Umklapp scattering is given by 

/2= -y|| r td'xcFr (A.l) 



Jo 



X (cos ( \/d>'KK(j){xi,Ti) + "^^^^^ j (JOS ( V^TTK(f){x2,T2) + ^^^^^ 



j' ^cos {V^cl>{x^,T,) + 



We want to evaluate this contribution here only in the thermodynamic limit L — » oo 
where it splits into a bulk and a boundary contribution. In this limit the two-point 
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correlation function for OBC at zero temperature is given by 

{(l){x-i,Ti)(l){X2,T2)) 



(A.2) 



1 

47r 



(In [{xt - X2f + - T2)2] - In [{xi + x^f + v\Tr - r^f]) . 



The exponentials of the bosonic field (j) can then by obtained by using the cumulant 

thooroni 



exp = exp {—A-KK{<j){x, t)(j){x, t))) 



and 



exp ^i\/87rK^^(xi, Ti)^ exp ^±iV87rif0(a;2, T2)^ 

= exp (-47rif [(02(a;i,Ti)) + (</.'(a;2, T2)) ± 2(</)(a;i,Ti)0(a;2,r2))]) 
Using the fact that (sin(V87ri4r^)) = 0, the susceptibility is given by 



(A.3) 
(A.4) 
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(A.5) 
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dx ( cos 



Using (A.3, A.4) the correlation functions can be evaluated. The usual conformal 
mapping then allows to obtain the second order correction at finite temperature 

V2 ir2 /^rp^ 4K-5 i-T^TL/v 

(A.6) 
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2K\ 



2K 



- 2{xi+X2f 



sinh(2xi) sinh(2x2) 



If wc arc far from the boundary, xi, X2 3> 1, then only the first term contributes and 
the integral gives a contribution to the bulk susceptibility 

^0 



L 

X (Xi - X2)^ 



d'x 



dr 



sinh(xi — X2 + ir) sinh(xi — X2 — ir) 



2K 



(A.7) 
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[i(cosh2x-cos2t)] 
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This integral can be evaluated analytically yielding the known result [4, 25] 
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To obtain the boundary contribution X2 have to subtract (A. 7) from (A. 6). As 
the kernel will then go to zero away from the boundary, we can shift the upper bound 
in the integration to infinity leading to 

4K- 



X2 



X2 



' L u3 



V 



dT 



X [sinh 2x sinh 2?/ (cosh2(x + y) — cos2t) (cosh2(a; — y) — cos2t)] 



-2K 



X {{xl+ xl) 



(cosh2(a; + y) — cos2t)^^ — (cosh2(a; — y) — cos2t) 



2K 



— [2 sinh 2x sinh 2y (cosh 2{x + y) — cos 2t)\ 



2K 



— 2xiX2 



(cosh2(a; + y) — cos2f)^^ — (cosh2(a; — y) — cos2f) 



4K 



— [2 sinh 2x sinh 2y (cosh 2{x + y) — cos 2t)]^^ | 



(A.9) 



The integral is convergent for ii' < 3/2 (A > 1/2) and can be evaluated numerically. 
Note, that this is the result for a semi- infinite chain, i.e., this result has to be multiplied 
by a factor 2 for an open chain segment. 



Appendix B. RG-improved coupling constants in the isotropic case 

The coupling constants in (5.1) fulfil the following set of RG equations [4] 
d9\\ 2gl dg± _ '^9\\9^ 



(B.l) 



rflnr 2 — 5f|| dlnr '^~9\\ 

Here r denotes the appropriate RG scale (T or L in the cases considered here). 
Following Lukyanov [4] the solution to this set of equations can be parametrised as 

5||=2(l-l/i^)l±^, 5^ = -4(l-l/i^)^ (B.2) 

where 

\2K-2 



,(1-,)— =^_j . (B.3) 

Here ro is a constant which can be chosen freely. 

This is the setup which quite generally allows to obtain results for the isotropic 
model by expanding the results for the anisotropic model in a power series in (l — l/K). 
From the scaling part (2.9) we obtain for the bulk susceptibility 

n=0 n=0 

Note, that the power series in produces terms higher order in q apart from the 
g°-term required. In particular, g\\/2 = (1 - 1/K){1 + 2q + 0{q^) with (1 - llK)q ~ 
(1 — 1/K)T'^^~^. Such a term is indeed present when expanding the second order 
contribution in A, Eq. (A. 8), in powers of (1 — l/K) and fixes the scale ro in (B.3). 
In other words, the scale in (B.3) is uniquely fixed because we have two equations for 
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the two unknowns: The pre-factor of g\\ is fixed through (B.4) and the expansion of 
(A. 8) fixes the scale to 

^ iK-i 



5(1 - q) 



2K-2 



re-V4^r(i + 2^) 



vT 1 



K 



2K-2 



r2(2 - K) 



(B.5) 



Comparing with the expression for the couphng constant A, Eq. (2.12), we see that 
we can write 

2K-2 

X=iK-l)^{l-qy 
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(B.6) 

(B.7) 



which in lowest order in (1 — 1/K) and q reduces to 

4 

Next, we consider the first order result for the boundary susceptibility (4.6). Using 
(B.6) we find 



6x1 = K^{1 - K)e^^-^^^^B{K, 1 - 2K) [n^ - 2i(;'{K)] 
Expanding this expression in (1 — 1/K) and q yields 



T5Xl = ^^^^^/9+(l-l/i^)^ 
For the first two terms we can write 

^ 12 8 



tp" 



(B.8) 



(B.9) 



(B.IO) 



however, this does not reproduce the third term in (B.IO) correctly. To obtain 
agreement also in order {1 — 1/ K)'^q^^^ the scale in (B.5) has to be changed. This leads 
to minimal improvements only, as this term is not only second order in (1 — 1/K) but 
also third order in y/q. We therefore keep the scale as determined from the expansion 
of the bulk susceptibility, Eq. (B.5), which is convenient to describe x{L,T) at the 
isotropic point. More important than the third term in (B.IO) is actually the second 
order contribution (A. 9) derived in the previous section. This contribution can be 
expressed as 

^ ,K^{K - l)2e^-i(7(l - g)2^-2(-2/) , (B.ll) 



where / denotes the integral in (A. 9) and the factor 2 has been added to account for 
the two boundaries. To lowest order in (1 — l/ZT) we therefore obtain 
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0.83. Finally, with g\\ 



(B.12) 
g and 



and a numerical cvahiation yields I{K = 1) f 
g±_ -g Eqs. (B.IO) and (B.12) yield (5.5). 

g as a function of temperature T can be found as follows: Using the first equation 
of (B.2) we can write q=[g-2{l- l/K)]/[g + 2{l - l/K)]. Now we can put this into 
q{l — q)^^~'^ and expand in (1 — 1/K) loading to 



g(i - q) 



2K-2 



1 + 



--+2 + 21n - 



21n(l - 1/K) 



{1-l/K). 



(B.13) 
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On the other hand, we can expand (B.5) leading to 



re-i/40Fr(i 



2K-2 



VT 1 



K 



2K-2 ^ 

« 1 + [1 + 21n2 + 21n(l - l/iiT) -47- 21n7r + 41n(2T/J)] . (B.14) 
Comparing (B.13) and (B.14) finally leads to (5.6). 
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